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Abstract:  Linear  dynamic  wave  and  diffusion  wave  to  arbitrarily  large  flow  increases.  As  wave  amplitude 

analytical  solutions  are  obtained  for  a  small,  abrupt  increases  the  monoclinal  rating  curve  diverges  from 

flow  increase  from  an  initial  to  a  higher  steady  flow.  that  for  a  linear  wave,  and  the  maximum  Froude  num- 

Equations  for  the  celerities  of  points  along  the  wave  ber  and  energy  gradient  along  the  profile  increase  and 

profiles  are  developed  from  the  solutions  and  related  move  toward  the  leading  edge.  A  monoclinal-diffu- 

to  the  kinematic  wave  and  dynamic  wave  celerities.  sion  solution  is  developed  for  the  diffusion  wave  equa- 

The  linear  solutions  are  compared  systematically  in  a  tions,  and  dynamic  wave-diffusion  wave  compari- 

series  of  case  studies  to  evaluate  the  differences  sons  are  made  over  a  range  of  amplitudes  with  the 

caused  by  inertia.  These  comparisons  use  the  celeri-  same  case  studies  used  for  linear  waves.  General 

ties  of  selected  profile  points,  the  paths  of  these  points  dimensionless  monoclinal-diffusion  profiles  exist  for 

on  the  x-t  plane,  and  complete  profiles  at  selected  each  depth  ratio  across  the  wave,  while  correspond- 

times,  indicating  general  agreement  between  the  solu-  ing  monoclinal  wave  profiles  exhibit  minor,  case- 

tions.  Initial  diffusion  wave  inaccuracies  persist  over  specific  Froude  number  dependence.  Inertial  effects  on 

relatively  short  time  and  distance  scales  that  increase  the  monoclinal  profiles  occur  near  the  leading  edge, 

with  both  the  wave  diffusion  coefficient  and  Froude  increase  with  the  wave  amplitude  and  Froude  number, 

number.  The  nonlinear  monoclinal  wave  solution  par-  and  are  responsible  for  the  differences  between  the 

allels  that  of  the  linear  dynamic  wave  but  is  applicable  dimensionless  profiles. 


Coven  Abrupt  river  wave  moving  downstream  in  the  Sacandaga  River  (New  York)  near  its 
confiuence  with  the  Hudson  River.  This  surge  was  formed  by  a  sudden  fiow  increase  of 
Stewarts  Bridge  Dam  severai  kiiometers  upstream.  (Photo  by  M.  Ferrick.) 


How  to  get  copies  of  CRREL  technical  publications: 

Department  of  Defense  personnel  and  contractors  may  order  reports  through  the  Defense  Technical  Information  Center; 
DTIC-BR  SUITE  0944 
8725  JOHN  J  KINGMAN  RD 
FT  BELVOIRVA  22060-6218 
Telephone  1  800  225  3842 

E-mail  help@dtic.mil 

msorders@dtic.mil 

WWW  http://www.dtic.dla.mil/ 

All  others  may  order  reports  through  the  National  Technical  Information  Sen/ice; 

NTIS 

5285  PORT  ROYAL  RD 
SPRINGFIELD  VA  22161 
Telephone  1  703  487  4650 

1  703  487  4639  (TDD  for  the  hearing-impaired) 

E-mail  orders@ntis.fedworld.gov 

WWW  http://www.fedworld.gov/ntis/ntishome.html 

A  complete  list  of  all  CRREL  technical  publications  is  available  from 
USACRREL  (CECRL-LP) 

72  LYME  RD 

HANOVER  NH  03755-1290 

Telephone  1  603  646  4338 

E-mail  techpubs@crrel.usace.army.mil 

For  information  on  all  aspects  of  the  Cold  Regions  Research  and  Engineering  Laboratory,  visit  our  World  Wide  Web  site: 
http://www.crrel.usace.army.mil 


CRREL  Report  98-1 


US  Army  Corps 
of  Engineers® 

Cold  Regions  Research  & 
Engineering  Laboratory 


Analysis  of  Linear  and  Monoclinal 
River  Wave  Solutions 

Michael  G.  Ferrick  and  Nicholas  J.  Goodman  January  1998 


Prepared  for 

OFFICE  OF  THE  CHIEF  OF  ENGINEERS 


Approved  for  public  release;  distribution  is  unlimited. 


PREFACE 


This  report  was  prepared  by  Michael  G.  Ferrick,  Hydrologist,  and  Nicholas  J. 
Goodman,  Physical  Science  Aid,  Geological  Sciences  Division,  Research  and  Engi¬ 
neering  Directorate,  U.S.  Army  Cold  Regions  Research  and  Engineering  Labora¬ 
tory,  Hanover,  New  Hampshire.  Funding  for  this  research  was  provided  by  DA 
Project  DT  08,  Bosnia  Support,  and  by  CWIS  Project  7VC103,  Water  Resources  of  Cold 
Regions. 

The  authors  thank  S.  Daly  and  Dr.  K.  O'Neill  of  CRREL,  K.  Newyear  of  the 
University  of  Washington,  Dr.  G.  Schohl  of  the  Tennessee  Valley  Authority,  and 
two  anonymous  reviewers  for  their  thoughtful  comments  that  helped  to  improve 
this  manuscript.  They  also  thank  Dr.  Y.  Nakano  of  CRREL  for  his  independent 
derivation  of  eq  19. 


CONTENTS 


Preface .  ii 

Nomenclature .  v 

Introduction .  1 

Linear  river  wave  equations .  2 

Solutions  for  linear  dynamic  and  linear  diffusion  waves .  4 

Nonlinear  monoclinal  and  monoclinal-diffusion  waves .  7 

Comparison  of  the  linear  solutions .  9 

Analysis  of  the  monoclinal  solutions .  14 

Conclusions .  23 

Literature  cited .  24 

Abstract .  27 

ILLUSTRATIONS 

Figure 

1.  Initial  and  upstream  boundary  conditions  of  the  linear  equations,  and 

dimensionless  (j)  profile  development  with  dimensional  distance  and  time .  4 

2.  Dimensionless  celerities  of  selected  linear  diffusion  wave  and  dynamic  wave 

profile  points .  10 

3.  Traces  on  the  x-t  plane  of  selected  linear  diffusion  wave  and  dynamic  wave 

profile  points .  11 

4.  Comparison  between  origin  positions  of  the  moving  coordinate  system  in 

the  linear  diffusion  wave  and  dynamic  wave  models  for  all  cases .  12 

5.  Linear  diffusion  wave  and  dynamic  wave  profiles  and  small-amplitude 

monoclinal  wave  profile  m .  13 

6.  Dimensionless  monoclinal  profile  celerity  and  overrun  discharge  as  a  func¬ 

tion  of  depth  ratio .  15 

7.  Evaluation  of  eq  41  as  a  function  of  depth  ratio  for  selected  dimensionless 

depths  on  the  monoclinal  profile .  16 

8.  Case  I  monoclinal  and  monoclinal-diffusion  profiles  at  two  distance  scales 

for  depth  ratios  ranging  between  1.1  and  50 .  16 

9.  Case  II  monoclinal  and  monoclinal-diffusion  profiles  at  two  distance  scales 

for  depth  ratios  ranging  between  1.1  and  50 .  17 

10.  Case  III  monoclinal  and  monoclinal-diffusion  profiles  at  two  distance  scales 

for  depth  ratios  ranging  between  1.1  and  5 .  17 

11.  Case  IV  monoclinal  and  monoclinal-diffusion  profiles  at  two  distance  scales 

for  depth  ratios  ranging  between  1.1  and  10 .  18 

12.  Case  V  monoclinal  and  monoclinal-diffusion  profiles  at  two  distance  scales 

for  depth  ratios  ranging  between  1.1  and  35 .  19 

13.  General  monoclinal-diffusion  profiles  at  two  dimensionless  distance  scales 

for  depth  ratios  ranging  between  1.1  and  50 .  19 

iii 


Figure 

14.  Monoclinal  profiles  for  all  cases  at  two  dimensionless  distance  scales  for 

depth  ratios  ranging  between  1.1  and  10 .  20 

15.  Dimensionless  monoclinal  and  steady  state  rating  curves  for  depth  ratios 

ranging  between  2  and  50 .  21 

16.  Energy  gradient  and  Froude  number  parameter  E  along  the  monoclinal 

wave  profile  for  selected  depth  ratios .  21 

17.  Energy  gradient  and  Froude  number  parameter  E  as  a  function  of  depth 

ratio  for  selected  dimensionless  depths  along  the  profile .  22 

18.  Dimensionless  depth  and  corresponding  maximum  E  as  functions  of  depth 

ratio .  22 


TABLE 

Table 

1.  Case  studies  used  to  compare  solutions .  9 


IV 


NOMENCLATURE 


B  monoclinal  wave  overrun  unit  discharge  (m^/s-m) 


c 

Cdif/  Cdyn 


Ck 

Co 


C+,C. 

c„c 

D 

E 

F,Fq 

F(x,t,^),  Fx,  ft 
G(x,t,^),  Gx,  Gt 
8 

hJi 
S(,  So 
t 

U 

V,  Vo,  Vf 


Vi 
Vt 
v',  y' 
xo, 


X 

X 

y,yo,yi 


y>^r  ^steady 


yi 

yr 

ycr 

Y 

<!),  (|)o/  <|)f 


X/  <t>xt/  <t>  XX 

A(x),  C(x), 
H,  z,  a,  P 

Ci,  Cl,  C2,  C3, 
C4,  C5,  Ce,  C7 


dimensionless  celerity 

celerity  of  a  point  on  the  linear  diffusion  wave,  dynamic  wave  solution 
profiles  (m/s) 

kinematic  wave  celerity  (m/s) 

celerity  of  a  disturbance  in  still  water  at  the  initial  uniform  flow  depth 
(m/s) 

dynamic  wave  celerity  in  the  downstream,  upstream  directions  (m/s) 
dimensionless,  dimensional  (m^^Vs)  Chezy  conveyance  coefficients 
noninertial  or  inertial  wave  diffusion  coefficient  (m^/ s) 
ratio  of  energy  gradient  to  bed  slope  gradient 
Froude  number,  initial  uniform  flow  Froude  number 
implicit  form  of  linear  dynamic  wave  solution  and  its  derivatives 
implicit  form  of  linear  diffusion  wave  solution  and  its  derivatives 
acceleration  due  to  gravity  (m/s) 

zeroth-order,  first-order  modified  Bessel  functions  of  the  first  kind 
energy  and  channel  bed  gradients 
time  (s) 

monoclinal  profile  celerity  (m/s) 

flow  velocity,  initial  uniform  velocity,  uniform  velocity  following  the 
wave  (m/s) 

small  increase  in  flow  velocity  from  the  initial  uniform  flow  (m/s) 

ratio  of  final  to  initial  uniform  flow  velocities 

velocity  derivative  (1/s),  depth  derivative  with  respect  to  X 

distance  scale,  dimensionless  distance  in  dimensionless  monoclinal 
wave  equation  (m) 

distance  along  the  chaimel  (m) 

distance  in  a  coordinate  system  moving  at  speed  U  with  the  monoclinal 
wave (m) 

flow  depth,  initial  uniform  depth,  and  uniform  depth  following  the 
wave (m) 

dimensionless  depth,  monoclinal  wave  unit  discharge,  and  steady  unit 
discharge 

small  increase  in  flow  depth  from  the  initial  uniform  flow  (m) 
ratio  of  final  to  initial  uniform  flow  depths 

depth  variable  representing  the  inertia  of  the  monoclinal  wave  (m) 
monoclinal  wave  depth  parameter  (m) 

dependent  variable,  initial  and  final  boundary  values  of  linear  wave 
equations 

first  and  second  derivatives  of  ^  with  respect  to  x  and  t 
parameter  of  linear  dynamic  wave  equation  (s) 
parameters  of  the  linear  d)mamic  wave  solution 

parameters  of  the  monoclinal  wave  solution 


V 


Analysis  of  Linear  and  Monoclinal  River  Wave  Solutions 


MICHAEL  G.  FERRICK  AND  NICHOLAS  J.  GOODMAN 

INTRODUCTION 

Much  effort  over  several  years  has  gone  into  the  development  and  application  of  numer¬ 
ical  models  for  imsteady  river  flow  problems.  In  particular,  one-dimensional  numerical 
models  that  solve  either  the  nonlinear  dynamic  wave  or  diffusion  wave  equations  have  been 
applied  over  wide  ranges  of  river  and  flow  conditions.  Criteria  for  choosing  a  model  from 
among  these  and  other  alternatives  for  given  conditions  have  been  obtained  from  analyses 
of  linearized  equations  (Ponce  and  Simons  1977,  Ponce  et  al.  1978,  Menendez  and  Norscini 
1982,  Kundzewicz  and  Dooge  1989).  Dimensionless  parameters  of  the  nonlinear  equations 
have  also  been  used  for  model  selection  (Woolhiser  and  Liggett  1967,  Ferrick  1985).  Several 
authors  have  treated  the  diffusion  wave-dynamic  wave  modeling  decision  as  a  choice  be¬ 
tween  simplicity  and  accuracy.  However,  Lighthill  and  Whitham  (1955)  argued  that 
dynamic  waves  are  subordinated  when  the  flow  is  "well  subcritical,"  making  the  character¬ 
istics  of  the  dynamic  wave  system  unsuitable  as  a  basis  for  computation.  Numerical 
dynamic  wave  and  diffusion  wave  models  cannot  be  readily  used  to  resolve  relative  accura¬ 
cy  issues  or  to  identify  optimal  model  selection  criteria. 

The  dynamic  wave  equations  include  flow  inertia  terms,  and  form  a  second-order  hyper¬ 
bolic  system  with  two  sets  of  characteristics  that  trace  the  paths  of  dynamic  waves  on  the 
x-t  plane.  The  diffusion  wave  equations  neglect  the  inertia  terms  as  small,  resulting  in  a  para¬ 
bolic  system  that  models  a  diffusing  "mass  wave."  As  the  magnitude  of  the  wave  diffusion 
term  decreases,  this  system  approaches  a  zero-diffusion  limit,  the  kinematic  wave  equation. 
This  first-order  hyperbolic  equation  has  a  single  set  of  characteristics,  the  subcharacteristics 
of  the  dynamic  wave  equations,  that  trace  the  paths  of  kinematic  waves  on  the  x-t  plane. 
Dynamic  waves  and  kinematic  waves  are  both  present  during  imsteady  river  flow,  and  it  is 
difficult  to  conceptualize  their  respective  roles  in  the  dynamic  wave  and  diffusion  wave 
models.  An  improved  understanding  of  these  models  would  be  an  important  step  toward 
resolution  of  relative  accuracy  and  model  selection  questions. 

Linearized  forms  of  the  dynamic  wave  and  diffusion  wave  equations  have  been  solved 
analytically  to  obtain  approximate  river  flows  (Dooge  and  Harley  1967,  Hayami  1951).  With 
variable  coefficients  treated  as  constants,  linear  solutions  are  strictly  valid  only  for  small 
flow  disturbances.  However,  these  solutions  are  valuable  because  of  their  common  structure 
with  corresponding  nonlinear  solutions.  The  linear  dynamic  wave  solution  is  the  most  gen¬ 
eral  and  provides  a  standard  for  comparison  with  simpler  linear  solutions,  but  systematic 
comparisons  have  not  been  developed.  Potential  benefits  include  better  definition  of  the 
correspondence  between  models,  and  resolution  of  the  time  and  distance  scales  where  dif¬ 
ferences  are  important.  Relationships  between  dynamic  wave  and  kinematic  wave  celerities 
and  the  downstream  translation  of  linear  wave  profiles  have  not  been  quantified  because 
equations  for  the  celerity  of  points  along  these  profiles  are  not  available.  Such  celerity  rela¬ 
tions  would  clarify  the  roles  of  characteristics  and  subcharacteristics  and  provide  insight 
into  the  structure  of  each  solution. 


Linear  dynamic  wave  and  diffusion  wave  solution  comparisons  cannot  fully  resolve  the 
relationship  between  these  models  because  large-amplitude  flow  increases  of  practical 
interest  must  be  described  by  the  nonlinear  equations.  This  deficiency  can  be  remedied  in 
part  by  considering  the  monoclinal  rising  wave,  a  nonlinear  dynamic  wave  analytical  solu¬ 
tion  (Chow  1959,  Henderson  1966,  Whitham  1974,  Hunt  1987,  and  Agsorn  and  Dooge 
1991).  The  monoclinal  wave  profile  is  an  arbitrarily  large  transition  between  low  steady, 
uniform  flow  downstream  and  high  steady  flow  upstream.  This  profile  represents  the  bal¬ 
ance  between  nonlinear  wave  steepening  and  diffusion,  and  has  a  known  constant  celerity 
that  increases  with  wave  amplitude  and  the  kinematic  wave  celerity.  A  comparison  of  a 
corresponding  nonlinear  diffusion  wave  equation  solution  with  the  monoclinal  wave 
would  identify  temporally  persistent  inertial  effects.  The  monoclinal  wave  solution  does 
not  describe  profile  development  nor  provide  the  time  and  travel  distance  required  to 
attain  a  steady  form.  However,  the  combination  of  linear  wave  and  monoclinal  wave  analy¬ 
ses  would  quantify  most  aspects  of  relative  dynamic  wave-diffusion  wave  solution 
behavior. 

The  purpose  of  this  report  is  to  utilize  analytical  solutions  to  better  understand  the  struc¬ 
ture  and  relative  behavior  of  the  dynamic  wave  and  diffusion  wave  unsteady  river  flow 
models.  An  abrupt  flow  increase  between  initial  and  final  steady  flows  is  used  as  an  up¬ 
stream  boundary  condition  to  maximize  the  contribution  of  inertia.  We  compare  linear 
dynamic  wave  and  diffusion  wave  solutions  in  a  series  of  subcritical  flow  case  studies. 
Equations  for  the  celerity  of  points  along  each  profile  are  derived  for  comparison  and  to 
explore  the  relationships  between  these  profile  celerities  and  the  dynamic  wave  and  kine¬ 
matic  wave  celerities.  This  development  provides  the  capability  to  trace  selected  profile 
points  on  the  x-t  plane,  another  means  to  compare  the  solutions.  We  also  compare  linear 
wave  profiles  to  depict  relative  behavior  through  time,  and  give  small-amplitude  mono¬ 
clinal  profiles  to  assess  progress  toward  equilibrium.  Our  nonlinear  monoclinal  wave  analysis 
uses  the  same  case  studies  as  for  linear  waves,  but  considers  a  range  of  wave  amplitudes.  The 
nonlinear  diffusion  wave  equations  are  solved  to  obtain  monoclinal-diffusion  profiles 
for  comparison  with  the  monoclinal  wave.  Relative  steepening  near  the  leading  edge  of  the 
monoclinal  profile  is  caused  by  flow  inertia  that  persists  through  time.  Nonlinear  effects 
increase  with  wave  amplitude,  progressively  separating  monoclinal  profile  shapes,  celeri¬ 
ties,  rating  curves,  and  the  Froude  numbers  and  flow  energy  gradients  along  the  profile 
from  those  of  linear  waves.  Dimensionless  dependent  variables  of  the  linear  and  monocli¬ 
nal  solutions  provide  ease  of  comparison,  while  dimensional  independent  variables  dis¬ 
tance  and  time  complement  physical  intuition.  For  generality  we  also  develop  and  com¬ 
pare  fully  dimensionless  monoclinal  and  monoclinal-diffusion  profiles. 


LINEAR  RIVER  WAVE  EQUATIONS 

The  continuity  and  momentum  equations  of  unsteady  flow  in  a  wide  rectangular  open 
channel  with  no  lateral  inflow  or  outflow  are  well  known  (Stoker  1957,  Mahmood  and 
Yevjevich  1975): 


dy  dy  dv  „ 

dt  dx  ^  dx 


(1) 


dv 

¥ 


dv  dy 


+  giS(-S^)  =  Q 


(2) 


where  depth  y  and  cross-sectional  average  velocity  v  are  dependent  variables,  Sf  is  the  flow 
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energy  gradient  obtained  from  the  Chezy  equation  with  a  constant  conveyance  coefficient, 
g  is  acceleration  due  to  gravity,  Sq  is  the  channel  bed  slope,  x  is  distance  and  t  is  time. 

Lighthill  and  Whitham  (1955)  obtained  a  linear  form  of  eq  1  and  2  by  substituting  for  y 
and  V 


y  =  yo+y\ 

v  =  Vq+v^  (3) 

where  constant  yo  and  Vq  represent  steady  uniform  flow  in  a  charmel  with  constant  slope 
and  resistance,  and  yj  and  Uj  represent  small  departures  from  that  flow  as 
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The  momentum  equation  (eq  5)  describes  linear  dynamic  waves.  If  we  assume  that  the 
inertia  terms  of  eq  5  are  small  relative  to  the  other  terms,  the  momentum  equation  for  linear 
diffusion  waves  is  obtained: 


3yi  ,  c:  f 
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(6) 


The  momentum  equation  for  linear  kinematic  waves  can  be  obtained  from  eq  6  by  assum¬ 
ing  the  depth  gradient  is  small  relative  to  the  bed  slope  and  can  be  neglected. 

Following  Lighthill  and  Whitham  (1955)  we  combine  eq  4  and  5  to  eliminate  either  the 
depth  or  velocity  derivatives.  The  resulting  second-order  linear  equations  for  depth  or  vel¬ 
ocity  are  the  same,  and  equivalent  to  the  original  system  of  first-order  equations: 


<l>tt  +  2^0 (t)xt  +(^0  ~^o)(^0  +Co)<t)xx +  “(<!> t  +  ^k'i’x)  -0 


(7) 


where  the  dependent  variable  (])  represents  either  Uj  or  yj,  subscripts  x  and  t  indicate  differ¬ 
entiation  with  respect  to  those  variables,  Cp  =  is  the  celerity  of  a  disturbance  in  still 
water,  =  3vo/2  is  the  kinematic  wave  celerity,  and  ti  =  Vo/2gSo.  Mendoza  (1995)  provided 
an  historical  perspective  on  the  development  and  solution  of  this  hyperbolic  equation, 
which  we  call  the  linear  dynamic  wave  equation.  Higher-order  dynamic  waves  travel 
along  two  sets  of  characteristics  described  by 


juo  +  co  =c+ 
dt  \uo  -  Co  =  c_ 


(8) 


where  the  product  of  c_  and  c+,  the  dynamic  wave  celerities  in  the  upstream  (-)  and  down¬ 
stream  (+)  directions  for  subcritical  flow,  appears  as  a  coefficient  in  eq  7. 

If  the  lower-order  terms  of  eq  7  were  absent,  corresponding  to  large  t],  the  general  solu¬ 
tion  would  have  the  form 

(t)  =  (!)i(A:-c+f)-i-(l)2(A:-c_f)  (9) 

with  a  structure  totally  dependent  on  the  dynamic  waves.  On  the  other  hand,  if  the  higher- 
order  terms  were  absent,  corresponding  to  r|  — >  0,  the  general  solution  would  have  the  form 

^  =  (foix-Cy,t)  (10) 
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depending  only  on  lower-order  kinematic  waves,  the  subcharacteristics  of  eq  7.  More  gen¬ 
erally  both  these  wave  types  occur  together,  with  small  disturbances  traveling  along  the 
characteristics  and  the  primary  mass  of  the  flow  moving  along  the  subcharacteristics. 

The  same  procedure  used  in  the  development  of  eq  7  can  be  repeated  with  eq  6  replacing 
eq  5,  yielding 

<t)t+C,,(l)^=D(j)xx  (11) 

where  D  =  CoTj  =  v^yQ  /  2So  is  a  noninertial  diffusion  coefficient.  We  call  this  parabolic  ad- 
vective-diffusion  equation  the  linear  diffusion  wave  equation.  The  dynamic  wave  celerities 
that  appear  explicitly  in  eq  7  depend  on  flow  inertia,  and  are  absent  from  eq  11.  The 
Hayami  (1951)  solution  of  eq  11  is  presented  in  Henderson  (1966).  Ferrick  et  al.  (1984)  de¬ 
veloped  an  equation  from  eq  7  with  the  same  form  as  eq  11  that  included  higher-order  x- 
derivatives.  The  presence  of  the  inertia  terms  resulted  in  an  inertial  diffusion  coefficient 
that  depends  on  the  Froude  number  Fq  -  vq/cq: 

D  =  c2il(l-Fo2/4).  (12) 

Dooge  (1973),  Whitham  (1974)  and  Menendez  (1993)  also  obtained  this  inertial  diffusion 
coefficient  using  different  methods. 


SOLUTIONS  FOR  LINEAR 
DYNAMIC  AND  LINEAR  ^  ° 
DIFFUSION  WAVES  o  ^ 

We  will  examine  the 
linear  dynamic  wave  and 
diffusion  wave  solutions 
for  an  instantaneous  in¬ 
crease  between  steady 
states  (t)o  and  ([)£  at  the  up¬ 
stream  boundary  {x  =  0) 
at  f  =  0.  Before  this  in¬ 
crease,  steady-uniform 
flow  conditions  down¬ 
stream  at  (t)o  are  assumed. 

This  initial  shock  maxi¬ 
mizes  the  importance  of 
inertia,  the  difference  be¬ 


X(km) 

Figure  1.  Initial  and  upstream  boundary  conditions  of  the  linear  equa¬ 
tions,  and  dimensionless  ^  profile  development  with  dimensional  dis¬ 
tance  and  time. 


tween  the  dynamic  wave  and  diffusion  wave  models.  The  solutions  ^{x,t)  are  required  for 
positive  X  and  t,  and  for  simplicity  we  define  dimensionless  ^  as 


$  = 


<l>-<|)o 
<l>f  -<1>0 


(13) 


with  $  =  1  at  the  upstream  boundary  and  a  solution  interval  of  zero  to  one.  In  the  remain¬ 
der  of  this  report  we  drop  the  tilde,  but  dimensionless  (])  is  implied.  Dimensional  {x,t)  and 
dimensionless  (])  are  used  in  Figure  1  to  depict  the  initial  and  upstream  bormdary  condi¬ 
tions,  and  profile  and  shock  development.  Lines  that  trace  the  motion  of  constant  ^  profile 
points,  and  profile  celerity  differences  between  these  points,  related  to  diffusion  of  the  pro¬ 
file,  are  indicated.  Nonlinear  wave  steepening  that  opposes  diffusion  can  cause  a  shock  to 
be  sustained,  but  it  is  not  present  here.  As  a  result,  the  imposed  shock  attenuates  and  the 
profile  separating  the  steady  states  elongates  with  time  and  distance. 
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The  Lighthill  and  Whitham  (1955)  subcritical  flow  (Fq  <  1)  solution  of  the  linear  dynamic 
wave  equation  (eq  7)  for  steady-uniform  initial  flow  with  a  unit-step  increase  in  ^  at  the 
upstream  boundary  is 


t-— 


0,f- 


where 


(14) 


a  = 


P  = 


(l  +  foV2) 


2q 


(1-F2)  (I-F0V4) 

2q 


1/2 


A{x)  = 


P  gfo^/4nco 

(l-foK 


1  fl-Fo/2^ 
2qcot  1  +  fo  y 


Ii  is  a  first-order  modified  Bessel  function  of  the  first  kind,  and  (j)(0,t)  is  the  upstream 
boimdary  condition  that  is  equal  to  0  for 


t- 


c+ 


<0 


and  equal  to  1  otherwise. 

Given  x  and  t  we  can  obtain  (j)  with  eq  14,  but  different  dependent  variable  designations 
would  be  helpful  for  studying  the  solution.  For  example,  insights  could  be  obtained  from 
wave  profiles  at  specified  times  and  from  constant  (j)  trajectory  traces  on  the  x-t  plane.  Also, 
it  would  be  useful  to  calculate  and  relate  the  celerity  of  a  point  on  the  dynamic  wave  profile 
to  the  dynamic  wave  and  kinematic  wave  celerities.  To  obtain  these  results  we  write  eq  14 
as  an  implicit  function: 


F{x,t,(Sf)  =  x-A{x)l  +  =  (15) 

0  ^  \  ‘^+J 


The  dynamic  wave  profile  celerity  of  a  point  with  constant  (j)  can  be  obtained  by  differenti¬ 
ating  F(x,f,(j))  and  setting  the  result  to  zero  as 


_dx  _  -Ff 


(16) 


where  Fx  and  Ft  are  the  partial  derivatives  of  F(x,f,^)  with  respect  to  x  and  f,  respectively. 
Dividing  eq  16  by  the  kinematic  wave  celerity  yields  dimensionless  profile  celerity 


^dyn 


^dyn 

Ck 


(17) 
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The  partial  derivative  in  eq  16  can  be  obtained  from  eq  15  as 


1+ 


LV 


FqX 

4tico 


fi(PzV^)  xdzf^ 

zi/2  za:rU  ''  J. 


e-ait-n  H 


where 


(18) 


and  Iq  is  a  zeroth-order  modified  Bessel  function  of  the  first  kind.  Similarly,  the  partial 
derivative  Ft  in  eq  16  can  be  obtained  from  eq  15  as 


Ft=x-A{x) 


CM+f 


ti(|}zi/^)  18z 

zl/2  Zdt 


Ja(Pz^Y 

zl/2  J. 


.  (19) 


In  eq  18  and  19 


X 

f-  — >0 

c+ 


and 


0,t-- 

V  C+y 


=  1 


are  assumed. 

F{x,t,^),  Fx  and  Ft  can  be  computed  using  subroutines  given  by  Press  et  al.  (1992)  for  10- 
point  Gauss-Legendre  integration  and  for  polynomial  approximation  of  modified  Bessel 
functions.  We  obtain  the  dynamic  wave  profile  at  selected  times  by  fixing  (])  and  f  and  calcu¬ 
lating  the  corresponding  x  using  eq  15  and  18  in  Newton's  method.  The  half-interval 
method  can  also  be  used  to  obtain  these  profiles  if  Newton's  method  fails  to  converge. 
Constant  ^  trajectory  traces  on  the  x-t  plane  are  obtained  by  specifying  (j)  and  x,  and  finding 
the  corresponding  t  using  Newton's  method  with  eq  15  and  19. 

Carslaw  and  Jaeger  (1959)  gave  the  solution  of  the  linear  diffusion  wave  equation  (eq  11) 
subject  to  the  initial  and  boimdary  conditions  stated  above,  and  we  write  it  as  an  implicit 
function: 


G(x,f,(t))  =  -(erfc| 


X-Ci-t 

(4Dt)i/2 


-exp  erfc 
D  J 


X  +  Ci^t 

(4D0^ 


-(t)(x,f)  =  0. 


(20) 


In  parallel  with  the  dynamic  wave  development  above,  for  a  given  constant  (j),  the  partial 
derivative  of  G{x,t,^)  with  respect  to  x  can  be  obtained  from  eq  20,  and  following  simplifica¬ 
tion  we  obtain 


^  -1 

Gx=7r;^7iT7Yexp| 


(KDf)i/2 


4Dt 


-^expf^l 

2D  ^  V  D  J 


erfc 


X  +  Ci^t 

(4Dfji^ 


The  partial  derivative  Gt  can  be  obtained  from  eq  20  and  simplified  as 

X 


Gt  = 


2t{nDt) 


^-(x-ci,tf^ 


V 


4Dt 


(21) 


(22) 
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We  use  the  same  development  as  above  to  obtain  the  dimensionless  profile  celerity  for  a 
fixed  value  of  (j)  on  the  diffusion  wave  profile,  which  in  simplified  form  is 


^dif 


-Gt  /  Gx 

Ck 


Ckf 


(23) 


The  forms  of  the  linear  diffusion  wave  and  dynamic  wave  solutions  have  little  apparent 
resemblance  to  each  other.  The  kinematic  wave  celerity  appears  in  each  term  of  the  diffu¬ 
sion  wave  solution  (eq  20),  but  not  in  the  dynamic  wave  solution  (eq  14).  The  dynamic 
wave  celerities  in  the  downstream  and  upstream  c_  directions  both  appear  in  z  of  the 
dynamic  wave  solution,  and  c+  provides  an  upper  bound  on  the  speed  of  disturbances 
moving  downstream.  There  is  no  similar  restriction  on  disturbances  traveling  downstream 
in  the  diffusion  wave  solution.  The  diffusion  coefficient  D  is  an  important  parameter  of  the 
diffusion  wave  solution,  while  rjCo  and  t|  are  corresponding  parameters  of  the  dynamic 
wave  solution.  The  Froude  number  Fq  appears  often  in  the  dynamic  wave  solution,  but  is 
absent  from  the  diffusion  wave  solution  imless  the  inertial  diffusion  coefficient  is  used. 


NONLINEAR  MONOCLINAL  AND  MONOCLINAL-DIFFUSION  WAVES 

Nonlinear  monoclinal  wave  solutions  that  are  analogous  to  the  linear  dynamic  wave 
and  diffusion  wave  solutions  will  now  be  developed  and  compared.  The  term  "monoclinal 
wave"  refers  to  the  classical  solution,  and  "monoclinal-diffusion  wave"  is  the  solution  de¬ 
veloped  after  neglecting  the  inertia  terms  of  momentum  equation  (eq  2).  A  monoclinal 
wave  profile  does  not  exist  for  the  kinematic  wave  equation  because  diffusion  is  not 
present  to  balance  nonlinear  steepening.  For  completeness  we  summarize  the  monoclinal 
wave  development  of  Whitham  (1974),  emphasizing  the  contribution  of  inertia  to  the  solu¬ 
tion. 

We  seek  a  solution  that  depends  on  a  single  variable  X  =  x-Ut,  where  U  is  the  constant 
profile  celerity,  and  rewrite  eq  1  and  2  as 

-Uy'  +  vy'  +  yv'  =  [y(LZ  -  u)]'  =  0  (24) 

{v-U)v'  +  gy'  =  gSQ--^  (25) 

Cty 

respectively,  where  Chezy  conveyance  coefficients,  dimensionless  and  dimensional  C, 
are  related  as 


The  first  term  on  the  left  side  of  eq  25  follows  from  the  inertia  terms  of  eq  2.  Integration  of 
eq  24  yields 

y(U-v)  =  B  (26) 

where  the  constant  of  integration  B  represents  a  wave  overrrm  unit  discharge. 

The  profile  celerity  U  is  obtained  from  eq  26  using  the  flow  states  on  either  side  of  the 
wave  as 


U  = 


ML1Mo.  =  („c?S  ]  "^0 

Vf-yo  °  V  yf -j/o  ^ 


=  Vi  + 


Vf+Vo 


(27) 
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where  the  alternate  forms  are  developed  by  eliminating  either  y  ox  v  using  the  Chezy  equa¬ 
tion.  Lf  in  eq  27  is  greater  than  Vi,  and  hence  all  velocities  along  the  profile.  With  eq  27  we 
rewrite  B  in  eq  26  as 


B  =  yoVi 


.yf-yo 


=  (^C2So) 


1/2 


ytyp 


=  igC}So)~' 


(  v}Vr  ^ 


1^0 


Vq+Vi 


Both  U  and  B  can  be  expressed  in  dimensionless  form  as 


(28) 


^^2 
Ck  3 


yr-1 


(29) 


B  ^yr(yr^^-l) 

mo  (yr-i) 


(30) 


These  parameters  increase  continuously  with  y^  from  minimums  of  1  and  1/2  at  =  1, 
respectively,  and  are  related  to  each  other  by  a  change  of  scale  as 


- 

Ck  3 


^oyo  j 


(31) 


In  eq  29  and  30  v,.  =  V(  /  Vq  can  be  substituted  for  y^  by  using  the  Chezy  equation  to  obtain 

(32) 

Returning  to  the  monoclinal  wave  equation  development,  we  divide  eq  25  by  gSo,  elim¬ 
inate  V  and  v'  from  the  inertia  term  using  eq  24  and  26,  and  obtain 


i: 

So 


B 


2  ^ 


(v/vof  . 
(y/yo) ' 


(33) 


Equation  26  is  used  again  to  eliminate  v  from  eq  33,  which  after  rearranging  becomes 

So 


_So(y-yf)(y-yo)(y-y) 


y^-ycr 


y  -ycr 


(34) 


where 


Y: 


yoyf 


<!/a 


and  y3.  =  B^/g  represents  the  contribution  of  inertia. 

Rewriting  eq  34  in  terms  of  convenient  dimensionless  variables  y*  =  y/yo  and  ^  =  x/xg 
yields 


dy* 

dx 


I  yo  J 


^  Y 

y*  — 
V  yo 


- 

y* 


.yo. 


(35) 
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where 


y  yt 
yo 


Vcc  _ 

yo 


fbyr 

yJ^^+1 


V/3 


Lighthill  and  Whitham  (1955)  deduced  the  monoclinal  profile  length  as  the  order  of  yo/So- 
Selecting  the  distance  scale  Xq  =  yo/So  simplifies  eq  35,  which  then  describes  dimensionless 
monoclinal  profiles  that  vary  with  y^  and  Fq,  and  dimensionless  monoclinal-diffusion  pro¬ 
files  that  depend  only  on  y,. 

The  solution  of  eq  34  for  the  monoclinal  wave  and  monoclinal-diffusion  wave  profiles 
can  be  obtained  by  separation  of  variables  and  integration  using  partial  fractions  as 

QX  +  Q  =  Cj  log(yf  -  y)  +  C3  log(y  -  yo)  +  C4  log(y  -  V) 


+Q(y{  -y)+Q(y“yo)+C7(y-Y) 


(36) 


where 

Q  =  So(yf  -  yo)  (yf  -  Y)  (yo  - ''") 
C2=(yf +ycr)(yo-^) 
C3=-(yo-ycr)(yf-‘'^) 
C4=(‘'"^-ya)(yf-yo) 

^5  =  ~yf  (yo“^) 
C6=-yo'(yf-v) 

Q  =^^(yf  “yo)- 


Cl  is  a  constant  that  results  from  the  integration  and  several  algebraic  manipulations.  We 
obtained  Ci  by  specifying  X  =  0  at  y  =  0.5(yf  +  yo).  Other  values  of  X  are  then  obtained  from 
eq  36  by  specifying  the  corresponding  y.  The  contribution  of  inertia  in  eq  36  is  eliminated  by 
setting  ycr  =  0,  and  the  monoclinal-diffusion  solution  results.  The  solution  for  v  is  obtained 
with  y  and  eq  26,  indicating  that  both  monoclinal  wave  types  have  the  same  rating  curve. 


COMPARISON  OF  THE 
LINEAR  SOLUTIONS 

The  five  cases  used  to  com¬ 
pare  the  linear  dynamic  wave 
and  diffusion  wave  solutions 
are  listed  in  Table  1.  These  cases 
represent  a  wide  range  of  con¬ 
ditions,  with  initial  velocities. 


Table  1.  Case  studies  used  to  compare  solutions. 

Vo  yo  T1  D  Co  Ck  c+ 


Case  (m/s) 

(m) 

So 

(s) 

(m2/s) 

(m/s) 

(m/s) 

(m/s) 

Fo 

I 

0.5 

1.0 

0.0005 

51 

500 

3.13 

0.75 

3.63 

0.16 

II 

0.5 

1.0 

0.0001 

255 

2500 

3.13 

0.75 

3.63 

0.16 

III 

3.0 

1.0 

0.0015 

102 

1000 

3.13 

4.5 

6.13 

0.96 

IV 

3.0 

3.0 

0.0005 

306 

9000 

5.42 

4.5 

8.42 

0.55 

V 

1.0 

3.0 

0.0001 

510 

15000 

5.42 

1.5 

6.42 

0.18 

initial  depths,  and  channel  slopes  that  vary  by  factors  of  6, 3  and  15,  respectively.  The  corre¬ 
sponding  variations  in  the  parameters  tj,  D,  and  Fq  by  factors  of  10, 30,  and  6,  are  also  large 
and  include  most  of  the  subcritical  flow  range.  Cases  I  and  II  represent  shallow,  low-veloc¬ 
ity  flows.  Parameters  ri,  D  and  Fq  each  have  relatively  small  values  in  case  I.  The  channel 
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slope  for  case  II  is  decreased  relative  to  case  I,  while  other  parameters  remain  unchanged, 
corresponding  to  significantly  increased  T]  and  D,  with  constant  Fq.  Case  III  is  a  shallow 
high  velocity  flow  with  channel  slope  and  velocity  increased  relative  to  the  first  two  cases. 
Parameters  r)  and  D  are  near  those  of  the  first  case,  but  Fq  is  at  the  high  end  of  the  range. 
Case  IV  represents  a  deeper  and  higher  velocity  flow  than  case  I  at  an  equal  channel  slope, 
corresponding  to  large  rj  and  D,  with  a  midrange  Fq.  Finally,  the  flow  in  case  V  is  deep  and 
relatively  slow  with  the  highest  ri  and  D,  and  Fq  near  that  of  cases  I  and  II. 

Dimensionless  profile  celerities  of  selected  points  ((j)  =  0.1, 0.3, 0.5, 0.7, 0.9)  were  obtained 
with  eq  17  and  23  for  dynamic  waves  and  diffusion  waves,  respectively,  and  results  for  all 
cases  are  presented  in  Figure  2.  In  early  time,  the  leading  edge  of  the  diffusion  wave  profile 
moves  downstream  at  a  celerity  that  initially  exceeds  and  later  is  less  than  the  dynamic 
wave  celerity  c+.  The  dynamic  wave  and  diffusion  wave  profile  celerities  of  all  points  con¬ 
verge  toward  each  other  and  the  kinematic  wave  celerity  with  time.  At  the  midpoint  (j)  =  0.5 
both  profile  celerities  rapidly  approach  the  kinematic  wave  celerity  and  then  remain  con¬ 
stant.  Smaller  ({)  profile  celerities  remain  higher  than  those  of  larger  (j),  and  diffusion  contin¬ 
ues  beyond  30  km  in  all  cases. 

In  case  I  the  discrepancies  between  diffusion  wave  and  dynamic  wave  profile  celerities 


Figure  2.  Dimensionless  celerities  of  selected  linear  diffusion  wave  and  dynamic  wave  profile  points 
for  all  cases  plus  case  III  with  the  modified  inertial  diffusion  coefficient.  Note  the  changes  in  the 
celerity  scale  between  panels. 
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are  minor  at  all  points  after  1  km  of  travel  distance.  Initial  differences  between  the  profile 
celerities  in  case  II  occur  as  a  result  of  the  shock  in  the  dynamic  wave  solution  for  (j)  up  to 
0.5,  and  some  differences  persist  for  5  km  downstream.  Similar  profile  celerity  disagree¬ 
ment  is  also  evident  in  case  III,  but  celerity  convergence  after  shock  attenuation  is  inexact, 
with  a  slightly  larger  range  remaining  in  the  diffusion  wave  solution.  With  the  modified 
inertial  diffusion  coefficient  (eq  12)  all  case  III  diffusion  wave  and  dynamic  wave  profile 
celerities  converge  by  6  km,  the  distance  for  shock  attenuation  below  (])  =  0.1.  The  large  D  in 
case  IV  causes  leading  edge  diffusion  wave  profile  celerities  to  greatly  exceed  the  dynamic 
wave  celerity.  The  shock  persists  for  17  km  downstream,  again  delaying  profile  celerity 
agreement.  The  profile  celerity  change  in  case  IV  with  the  inertial  diffusion  coefficient  was 
negligible.  The  celerity  comparisons  for  case  V  are  similar  to  those  of  case  IV,  with  larger 
initial  dimensionless  profile  celerities  in  both  solutions,  and  celerity  agreement  at  all  points 
following  shock  attenuation  at  17  km. 

Traces  on  the  x-t  plane  of  selected  wave  profile  points  in  Figure  3  help  to  visualize  the 
effects  of  profile  celerity  differences.  The  time  scales  used  for  each  case  are  related  by  the 
kinematic  wave  celerities  given  in  Table  1.  The  dynamic  wave  /from  the  origin  at  f  =  0, 
termed  the  dynamic  forerunner  by  Stoker  (1957),  carries  the  initial  shock  downstream  at  a 
constant  celerity  c+.  A  positive  value  of  F{x,t,^)  immediately  behind  the  forerunner  indi¬ 
cates  that  a  given  (j)  is  on  the  forerunner.  The  (t>-traces  that  successively  separate  from  the 


Figure  3.  Traces  on  the  x-t  plane  of  selected  linear  diffusion  wave  and  dynamic  wave  profile  points  and 
the  dynamic  forerunner  ffor  all  cases  plus  case  III  with  the  modified  inertial  diffusion  coefficient.  Note 
the  changes  in  the  time  scale  between  panels. 
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dynamic  forerunner  indicate  a  progressively  diminishing  shock  amplitude.  Afterward,/ sep¬ 
arates  from  the  profile  and  no  longer  contributes  to  the  solution.  Overall,  the  x-t  traces  indi¬ 
cate  general  agreement  between  the  dynamic  wave  and  diffusion  wave  solutions  following 
attenuation  of  the  initial  shock. 

The  case  I  traces  of  all  corresponding  dynamic  wave  and  diffusion  wave  profile  points  are 
essentially  identical.  The  dynamic  forerurmer/in  case  I  begins  to  separate  from  the  profile  in 
the  first  5  km,  and  leads  by  increasing  distances  farther  downstream.  At  early  time  in  case  II 
the  diffusion  profile  celerity  of  the  leading  edge  exceeds  that  of/,  causing  minor  differences 
between  the  traces.  These  differences  disappear  about  5  km  from  the  origin,  and  thereafter  the 
traces  of  all  points  are  identical.  The  forerunner  and  profile  in  case  II  progressively  separate 
beyond  10  km  from  the  origin.  The  high  Fq  in  case  III  reduces  the  rate  of  spread  of  the  dynamic 
profile,  causing  the  diffusion  traces  of  the  front  half  to  lead  and  of  the  back  half  to  lag  the 
d3mamic  traces,  and  these  trends  persist.  With  the  inertial  diffusion  coefficient  the  traces  in 
case  III  agree  closely  beyond  8  km,  and  the  profile  progressively  separates  from /beyond  15 
km.  In  case  IV  a  shock  amplitude  of  0.1  is  carried  by  the  forerunner  for  16  km,  with/separat- 
ing  from  the  profile  beyond  25  km.  Excess  profile  spreading  in  the  diffusion  solution,  caused 
by  initial  celerity  differences  and  high  diffusion  coefficient,  persists  throughout  and  was  not 
greatly  improved  by  the  inertial  diffusion  coefficient.  The  case  V  dynamic  wave  and  diffusion 
wave  traces  compare  similarly  to  those  of  case  IV,  except  that  differences  higher  on  the  profile 
do  not  persist. 

We  can  compare  dynamic  wave  and  diffusion  wave  profiles  through  time  on  a  single  figure 
by  using  a  moving  x-coordinate  system  with  origin  at  (j)  =  0.5.  Comparisons  of  these  origin  x- 
values  through  time  (x^yn,  ;Cdif)  ^re  given  in  Figure  4.  The  origin  traces  of  the  dynamic  wave 
model  were  similar  for  cases  with  the  same  kinematic  wave  celerity.  The  absolute  value  of  the 
difference  in  origin  position  between  the  dynamic  wave  and  diffusion  wave  models  was  al- 


t(sx10-3) 


Figure  4.  Comparison  between  origin  positions  of  the  moving  coordinate  system  in  the  linear  diffusion 
wave  (xdip  and  dynamic  wave  (x^yn)  models  for  all  cases. 
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ways  less  than  800  m,  and  diminished  to  less  than  4  %  of  after  1  hour.  The  ratio  of  these 
origin  positions  rapidly  approached  1  from  below  in  all  cases.  Larger  diffusion  coefficients 
correspond  to  larger  absolute  differences  in  origin  position  early  and  generally  larger  ratio 
differences  from  1. 

Dynamic  wave  and  diffusion  wave  profiles  at  selected  times  are  compared  in  the  moving 
coordinate  system  in  Figure  5.  The  dynamic  wave  solution  includes  the  initial  shock  on  the 
forerunner.  At  early  time  the  leading  edge  of  the  diffusion  profile  always  precedes  that  of  the 
dynamic  profile  by  a  distance  that  increases  with  D  or  r)  and  the  Froude  number.  The  profiles 
in  case  I  rapidly  converge  and  are  nearly  identical  after  1200  s.  In  case  II  the  shock  front  is 
preserved  for  a  longer  time,  and  the  profiles  converge  by  10,800  s.  With  high  Fq  in  case  III  the 
profiles  tend  to  converge  after  the  shock  diminishes,  but  the  dynamic  wave  profile  retains 
more  steepness  than  the  diffusion  wave  profile.  Case  III  profiles  with  the  inertial  diffusion 
coefficient  are  nearly  identical  after  3600  s.  In  case  IV  the  shock  persists  for  a  longer  time,  and 
profile  convergence  requires  more  than  10,800  s.  The  minor  differences  in  profile  steepness 
remaining  at  large  times  can  be  minimized  with  the  inertial  diffusion  coefficient.  The  case  V 
profile  comparisons  are  similar  to  those  of  case  IV,  but  with  a  low  Froude  number  the  inertial 
diffusion  coefficient  is  not  needed  for  agreement  at  large  times.  Monoclinal  wave  profiles  of 


Figure  5.  Linear  dijfusion  wave  and  dynamic  wave  profiles  and  small-amplitude  monoclinal  wave  profile  mfor  all 
cases  plus  case  III  with  the  modified  inertial  dijfusion  coefficient.  Note  the  changes  in  the  distance  scale  between  panels. 
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amplitude  O.lyo,  representing  fully  diffused  linear  wave  profiles  at  large  times,  are  also 
given  for  each  case  in  Figure  5.  The  most  prominent  feature  of  these  solutions  is  the  low 
profile  slope,  and  the  extended  time  indicated  for  the  linear  wave  to  attain  this  profile.  Only 
the  higher  Froude  number  linear  profiles  at  36,000  s  even  approach  the  small-amplitude 
monoclinal  profiles. 


ANALYSIS  OF  THE  MONOCLINAL  SOLUTIONS 

The  numerator  of  the  depth  gradient  in  eq  34  does  not  change  sign  along  the  profile,  but 
a  sign  change  in  the  denominator  can  result  from  the  presence  of  inertia.  This  sign  change 
indicates  a  monoclinal  wave  profile  that  turns  back  upstream,  becomes  unstable,  and  forms 
a  shock.  In  contrast,  the  monoclinal-diffusion  profile  cannot  become  unstable  because  y^t  is 
not  present.  Initial  monoclinal  wave  instability  occurs  at  the  toe  of  the  profile  when  the 
denominator  goes  to  zero,  and 


y  =  yo  =  ycx 


s . 


(37) 


Following  Hunt  (1987)  we  evaluate  the  stability  limit  eq  37  using  eq  28,  and  after  some 
algebra,  the  depth  ratio  across  the  wave  is  obtained  as  a  function  of  Fq 


(38) 


The  depth  ratio  range  of  stable  monoclinal  wave  profiles  decreases  as  Froude  number  in¬ 
creases  toward  2.  Conversely  Fq  can  be  obtained  as  a  function  of  y^  with  eq  35  as 


(39) 


The  stability  limit  in  eq  37,  evaluated  using  eq  26,  yields  U  =  vq+cq  as  the  maximum 
profile  celerity  prior  to  instability.  The  range  of  stable  profile  celerities,  bounded  below  by 
the  celerity  of  lower-order  kinematic  waves  and  bounded  above  by  the  celerity  of  higher- 
order  dynamic  waves,  can  be  written  in  dimensionless  form  as 


2 

<- 

.3 


(40) 


Dimensionless  profile  celerity  and  overrun  discharge  given  in  Figure  6  increase  continu¬ 
ously  with  wave  amplitude  from  lower  limits  of  1  and  0.5,  respectively,  at  y^  =  ytlyo  =  1.  The 
upper  limits  are  indicated  for  selected  values  of  Fq  by  dots  that  follow  from  eq  40. 

The  difference  between  the  monoclinal  and  monoclinal-diffusion  solutions  can  be  para¬ 
meterized  by  considering  the  bracketed  term  in  eq  33  rewritten  using  dimensionless  vari¬ 
ables 


■  1- 


'yciif 

.  y . 


=  1- 


[(yr-i)y+if 


(41) 


where  dimensionless  depth  y,  defined  analogously  to  $  in  eq  13,  varies  between  0  and  1. 
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Figure  6.  Dimensionless  monoclinal  profile  celerity  and  overrun  discharge 
as  a  function  of  depth  ratio.  The  stable  profile  range  spans  a  depth  ratio  range 
from  1  to  an  upper  bound,  indicated  by  the  dots,  that  decreases  as  the  Froude 
number  increases.  Panel  b  is  an  expanded  view  of  the  shaded  area  of  panel  a. 


For  monoclinal-diffusion  waves  eq  41  has  a  value  of  1,  and  the  deviation  from  1  indicates 
the  relative  importance  of  inertia.  The  same  cases  analyzed  for  linear  waves  are  used  to 
depict  monoclinal  waves,  with  time  deleted  from  the  parameters  considered  and  depth 
ratio  across  the  wave  representing  amplitude  added.  Evaluations  of  eq  41  for  each  case  are 
presented  in  Figure  7  as  a  function  of  y,  for  selected  values  of  y ,  all  with  limit  1  -  Ff!  4  at  yr 
=  1.  At  low  Froude  numbers  the  part  of  the  profile  affected  by  inertia  is  very  close  to  leading 
edge  (small  y ),  and  then  only  when  depth  ratios  are  large.  Differences  between  monoclinal 
and  monoclinal-diffusion  profiles  near  the  leading  edge  increase  with  y^  and  Fq-  Negative 
values  of  eq  41  indicate  that  the  profile  point  y  is  located  on  a  shock.  Cases  I  and  II  have  a 
calculated  stability  limit  of  yr  =  51,  where  the  dimensionless  shock  amplitude  is  smaller 
than  0.001.  In  case  III  with  high  Fq  much  more  of  the  profile  is  affected  by  inertia,  and  larger 
shocks  occur  at  relatively  small  depth  ratios.  Case  IV  is  intermediate  between  these  condi¬ 
tions,  and  case  V  is  similar  to  cases  I  and  II. 

Monoclinal  and  monoclinal-diffusion  dimensionless  depth  profiles  for  case  I,  presented 
in  Figure  8,  are  in  exact  agreement  except  for  the  leading  edge  at  yr  =  50,  where  the  diffusion 
solution  leads.  The  front  half  of  these  profiles  shorten  and  steepen  as  yr  increases  to  10.  At 
yr  =  50  the  wave  front  lengthens,  and  the  steepest  portion  continues  forward  to  the  leading 
edge.  The  profile  comparisons  and  trends  for  case  II  in  Figure  9  are  identical  to  those  of  case 
I,  except  that  profile  lengths  are  significantly  increased  as  a  result  of  much  higher  diffusion. 
Case  III,  depicted  in  Figure  10,  has  monoclinal  and  monoclinal-diffusion  profiles  that  pro¬ 
gressively  separate  below  a  dimensionless  depth  of  about  0.4.  The  leading  edge  of  the 
steeper  monoclinal  profile  lags  behind  that  of  the  diffusion  profile.  At  y^  =  5,  outside  the 
stable  profile  range,  overrun  of  the  leading  edge  of  the  monoclinal  profile  indicates  shock 
formation  up  to  a  dimensionless  depth  of  about  0.1.  Case  IV  in  Figure  11  is  qualitatively 
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Figure  10.  Case  III  monoclinal  and  monoclinal-dijfusion  profiles  at 
two  distance  scales  for  depth  ratios  ranging  between  1.1  and  5. 
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Figure  11.  Case  IV  monoclinal  and  monoclinal-dijfusion  profiles  at 
two  distance  scales  for  depth  ratios  ranging  between  1.1  and  10. 


similar  to  case  III  with  distance  scales  substantially  increased.  At  yr  =  2  the  profiles  begin  to 
separate  at  a  dimensiorUess  depth  of  about  0.3,  and  an  overrun  of  the  leading  edge  occurs  at 
y,  =  10  up  to  a  dimensionless  depth  of  0.03.  Case  V,  presented  in  Figure  12,  is  qualitatively 
similar  to  cases  I  and  11  except  that  the  largest  diffusion  coefficient  produces  the  longest  pro¬ 
files  of  all  the  cases. 

General  results  from  these  comparisons  are  that  monoclinal  and  monoclinal-diffusion 
profiles  agree  for  all  values  of  y^  with  Fq  <  0.2,  and  that  profile  length  increases  with  T)  or  D. 
At  small  Fq  the  depth  ratio  needed  to  produce  a  shock  is  large,  and  the  shock  height  and 
overrun  distance  of  the  leading  edge  are  small.  These  results  agree  with  the  Lighthill  and 
Whitham  (1955)  contention  that  dynamic  waves  are  subordinated  at  "well  subcritical" 
Froude  numbers.  As  Fq  increases  the  monoclinal  waves  differ  over  a  larger  portion  of  the 
profile,  shocks  occur  at  smaller  depth  ratios  and  their  dimensionless  amplitudes  increase, 
and  for  a  given  T]  the  profile  length  decreases.  General  dimensionless  monoclinal-diffusion 
profiles  for  each  depth  ratio  are  given  in  Figure  13  as  a  function  of  x,  and  include  the  profiles 
of  all  cases  as  indicated  by  eq  35.  Similar  dimensionless  monoclinal  profile  plots  in  Figure  14 
are  almost  as  well-behaved,  but  differences  near  the  leading  edge  occur  for  each  y,.  due  to 
their  Fq  dependence. 

Steady  flow  rating  curves  relate  river  stage  or  mean  depth  at  a  given  location  to  a  unique 
discharge.  The  governing  equation  for  linear  waves  holds  with  either  u  or  y  as  the  depen¬ 
dent  variable,  indicating  a  fixed  steady  flow  rating.  In  unsteady  flow  the  discharge  relating 
to  a  given  stage  generally  varies  from  that  for  steady  flow,  depending  on  the  rate  of  rise  or 
fall  of  the  hydrograph.  We  will  develop  and  compare  dimensionless  monoclinal  wave  and 
steady  flow  ratings.  Using  eq  26,  an  equation  for  the  monoclinal  wave  unit  discharge  can  be 
written  in  terms  of  depth  and  depth  ratio  as 
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Figure  12.  Case  V  tnonoclinal  and  monoclinal-diffusion  profiles  at 
two  distance  scales  for  depth  ratios  ranging  between  1.1  and  35. 


Figure  13.  General  monoclinal- 
diffusion  profiles  at  two  dimen¬ 
sionless  distance  scales  for  depth 
ratios  ranging  between  1.1  and  50. 


X 


Figure  14.  Monoclinal  profiles  for  all  cases  at  two  dimensionless  distance  scales  for  depth 
ratios  ranging  between  1.1  and  10. 


(j  _  vy  _(y/yo)(yr^^-i)-yr(y?^^-i) 

%  Voyo  (yr-1) 


(42) 


Defining  dimensionless  imit  discharge  with  the  same  form  as  dimensionless  depth,  eq  2 
becomes 


q  =  y  (43) 

independent  of  the  depth  ratio.  Using  the  Chezy  equation  we  obtain  a  corresponding  rela¬ 
tionship  for  dimensionless  unit  discharge  in  steady  flow  conditions  as  a  function  of  depth 
ratio 


*7steady 


(44) 


In  the  linear  wave  limit  as  approaches  1,  the  rating  curves  represented  by  eq  43  and  44  are 
identical.  These  rating  curves,  given  in  Figure  15,  indicate  that  for  a  given  dimensionless 
depth  as  y^.  increases  the  monoclinal  wave  emit  discharge  also  increases  relative  to  that  for 
steady  flow. 

Froude  number  and  energy  gradient  along  a  linear  wave  are  unchanged  from  those  of 
the  initial  steady  flow.  Large  amplitude  monoclinal  waves  with  discharges  along  the  pro¬ 
file  that  greatly  exceed  those  of  steady  flow  at  comparable  depths  can  have  significantly 
larger  Sf  and  F.  We  define  E  =  Sf/Sg,  and  with  the  Chezy  equation  and  eq  26  obtain 


y/yo  UoJ  [(y,-l)y+lf  ■ 


(45) 
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Figure  15.  Dimensionless  monoclinal  and  steady  state  rating  curves  for 
depth  ratios  ranging  between  2  and  50. 


Figure  16  gives  E  as  a  function  of  y  for  selected  values  of  y^.  The  maximum  £  increases  with 
1/r  and  moves  from  near  the  midpoint  of  the  profile  toward  the  leading  edge.  Conversely, 
Figure  17  gives  £  as  a  function  of  y^  for  selected  values  of  y.  At  small  y^  the  largest  £  is 
located  near  the  midpoint  of  the  profile.  As  y^  increases  £  approaches  a  constant  starting  at 
the  midpoint  of  the  profile  and  progressing  toward  the  leading  edge.  At  large  depth  ratios 
£  is  inversely  proportional  to  g,  increasing  toward  the  leading  edge. 


Figure  16.  Energy  gradient  and  Froude  number  parameter  E  along  the  monocli¬ 
nal  wave  profile  for  selected  depth  ratios.  Panel  b  is  an  expanded  view  of  the  shaded 
area  in  panel  a. 
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The  maximum  E  is  then  obtained  by  substituting  eq  46  into  45  and  simplifying 

^  4  (yP-if 

The  location  and  value  of  En^^x  are  given  in  Figure  18.  E^gx  increases  nearly  linearly  over 
the  i/r  range,  and  its  location  rapidly  approaches  the  leading  edge  as  increases  from  1. 
Both  energy  gradient  and  Froude  number  vary  continuously  along  the  monoclinal  wave 
profile,  with  amplitudes  proportional  to  y^.  Froude  numbers  that  exceed  1  can  occur  on 
stable  profiles  when  Eq  is  large,  as  in  case  III. 

CONCLUSIONS 

The  presence  or  absence  of  the  inertia  terms  distinguishes  dynamic  wave  and  diffusion 
wave  models  of  unsteady  river  flow.  Analytical  solutions  of  the  linear  dynamic  wave  and 
diffusion  wave  equations  were  compared  for  a  small  instantaneous  increase  from  an  initial 
steady,  uniform  flow  condition  throughout  the  channel  to  a  higher  steady  flow  velocity, 
and  depth  at  the  upstream  boundary.  The  comparison  used  case  studies  that  represented  a 
wide  range  of  flow  depth,  velocity,  channel  slope,  and  wave  diffusion  coefficient,  and 
sparmed  the  range  of  subcritical  Froude  numbers.  Analytical  solutions  were  also  obtained 
for  nonlinear  monoclinal  wave  and  monoclinal-diffusion  wave  equations,  and  compari¬ 
sons  were  again  made  using  the  same  case  studies  with  a  wide  range  of  wave  amplitudes. 
The  linear  solution  comparisons  focused  on  the  evolution  of  the  dynamic  wave  and  diffu¬ 
sion  wave  profiles  with  time  and  distance,  while  the  nonlinear  solution  comparisons  inves¬ 
tigated  the  effects  of  wave  amplitude  and  persistent  inertia. 

The  linear  solution  comparisons  included  the  celerity  and  the  trace  on  the  x-t  plane  of 
selected  points  from  the  wave  profile,  and  the  complete  dynamic  and  diffusion  profiles  at 
selected  times.  The  initial  shock  traveled  downstream  with  the  dynamic  forerunner  at  c+, 
the  maximum  celerity  in  subcritical  flow.  A  limitation  of  the  diffusion  wave  solution  is 
premature  replacement  of  this  shock  by  a  profile  having  a  range  of  point  celerities  that 
exceed  c+  near  the  leading  edge.  The  diffusion  wave  and  dynamic  wave  profiles  remain 
distinct  imtil  after  the  shock  attenuates  and  their  profile  celerities  converge.  In  cases  where 
the  Froude  number  approaches  1,  this  convergence  requires  a  diffusion  coefficient  cor¬ 
rected  for  inertia.  Points  near  the  leading  edge  of  each  profile  travel  faster  than  those  higher 
on  the  profile,  causing  diffusion.  These  differences  diminish  over  time  and  distance,  and  all 
dynamic  and  diffusion  profile  celerities  asymptotically  approach  that  of  a  kinematic  wave. 
General  agreement  of  the  linear  diffusion  wave  and  dynamic  wave  solutions  after  attenua¬ 
tion  of  the  shock  is  indicated  in  all  cases  by  common  profile  celerities  and  x-t  traces,  and  by 
profile  covergence  with  time.  The  role  of  the  characteristics  in  the  linear  solution  becomes 
negligible  following  shock  attenuation  at  time  and  distance  scales  that  increase  with  both  r] 
and  the  Froude  number. 

The  analysis  of  the  nonlinear  monoclinal  wave  solutions  linked  important  inertial  effects 
at  large  time  with  increasing  Froude  number.  As  Fq  increases  corresponding  monoclinal 
and  monoclinal-diffusion  profiles  separate  near  the  leading  edge,  and  these  differences 
increase  with  wave  amplitude.  Monoclinal  profile  instability  occurs  at  higher  Eq,  but  mono¬ 
clinal-diffusion  profiles  are  always  stable.  General  dimensionless  monoclinal-diffusion 
profiles  exist  for  each  depth  ratio  with  distance  scaled  by  yo/So,  but  monoclinal  profiles 
deviate  from  these  general  profiles  at  higher  Eq.  Several  effects  of  wave  amplitude  on  mon¬ 
oclinal  waves  were  also  identified.  The  celerity  of  a  small-amplitude  monoclinal  wave 
equals  that  of  a  kinematic  wave,  and  it  increases  continuously  from  this  lower  limit  with 
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wave  amplitude.  The  steepness  near  the  front  of  the  monodinal  wave  profile  and  the  differ¬ 
ence  between  monodinal  and  steady  flow  or  linear  wave  rating  curves  both  increase  as 
amplitude  increases.  The  energy  gradient  and  Froude  number  at  all  points  along  the  profile 
also  increase  with  wave  amplitude,  and  the  location  of  the  maximum  shifts  continuously 
toward  the  leading  edge. 
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